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Abstract. Specific heat measurements have been successfully used to probe 
unconventional superconducting phases in one-band heavy-fermion and organic 
superconductors. We extend the method to study successive phase transitions 
in multi-band materials such as iron based superconductors. The signatures are 
multiple peaks in the specific heat, at low temperatures and high magnetic field, 
which can lead the experimental verification of unconventional superconducting 
states with non-zero total momentum. 
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1. Introduction 

Ideal diamagnetism is one of most striking properties of superconductors, which is 
manifested by the Meissner effect, the complete expulsion of a magnetic field from 
the volume of a superconductor. Conversely, magnetic fields with relatively large 
value can destroy superconductivity. In most real (II type) superconductors it may 
take place in two ways - by orbital or paramagnetic effects. The orbital pair¬ 
breaking is connected with the rise of the Abrikosov vortex state in superconductors, 
while paramagnetic pair-breaking originates from the Zeeman splitting of electronic 
energy levels. Both effects determine the upper critical magnetic field, in which the 
relative importance of the orbital and paramagnetic effects in the suppression of the 
superconductivity is described by the Maki parameter a = \/2H° 2 b / H^ 2 [1], the ratio 
of the critical magnetic fields at zero temperature H°'iff and H^ 2 , derived from orbital 
and diamagnetic effects respectively. 

In most superconductors, the orbital pair-breaking effects are more disruptive 
than the diamagnetic ones, thus H^ 2 is usually larger than H° 2 b (a <C 1) and 
superconductivity disappears when vortex cores begin to overlap. When a > 1, 
superconductivity is destroyed by the Zeeman effect (Pauli paramagnetism), and 
systems exhibiting this property are called Pauli limited materials. 

In absence of an external magnetic field, superconductors are found in the 
Bardeen-Cooper-Shrieffer (BCS) state [2, 3], where superconductivity is formed 
by Cooper pairs with total momentum equal zero. However, in Pauli-limited 
superconductors, the external field can lead to interesting phenomena (near or) above 
the critical magnetic field H^ 2 , such as the transition to the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase [4,5]. In contrast to the BCS state, now the Cooper 
pairs are formed between the spin-up and spin-down sheets of the split Fermi surface 
with non-zero total momentum. Moreover this phase exhibits a spatially oscillating 
superconducting order parameter (SOP) in real space and spin polarization. 

Although the FFLO phase was already theoretically described in the 1960s, 
the experimental search is still ongoing and riddled with difficulties. They come 
as a consequence of the physical properties of this inhomogeneous phase: it 
can occur only in Pauli limited superconductors at low temperature and high 
magnetic field (LTHM) regime. Only in the last decade systems have appeared in 
which we expect an experimental verification of this phase, such as heavy-fermion 
superconductors [6] (e.g. CeCoIn 5 [7-11]) or organic superconductors (e.g. /3"- 

(ET) 2 SF 5 CH 2 CF 2 S03 [12], A-(BETS) 2 FeCl 4 [13], k-(BEDT-TTFS) 2 Cu(NCS) 2 [14- 
18]). These chemical compounds are often modeled as effectively quasi-2T) one-band 
systems. However in 2013, D. Zocco et al. in Ref. [19] reported possibly the FFLO 
phase in the multi-band iron-based superconductor KFe 2 As 2 . Although it should 
be noted that hints about the FFLO phase in pnictides had been reported in other 
experimental [20-25] and theoretical [26] works previously. 

In experiments on the ordered phases (such as superconducting or magnetic), 
measurements of the anomalies in the specific heat C are one of the most sensitive 
tools [27]. Using this method we can study phase transitions and their type, or the 
nodes of the gap function in the superconducting state [28-30] . In this sense peaks in 
the specific heat reveal information about the phase transition. Moreover, a narrow 
peak in C(T) is associated with a first order transition, while the A-shape behavior 
is typical for the second order transition [31,32]. This method can be used to find 
phase transition e.g. in multiband superconductors like IBSC [33-36] or MgB 2 [37,38]. 
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Specific heat measurements have been also successfully used to investigate the FFLO 
phase in heavy-fermion systems [7,8,11] and organic superconductors [14], where the 
peaks, deep in the superconducting state in the LTHM regime, have been interpreted 
as phase transitions from BCS to FFLO state. 

Because the experimental evidence to confirm the FFLO phase in iron-based 
superconductors (IBSC) is still needed, the description of the theoretical properties of 
this phase is paramount. In this paper we discuss experimental consequences multiple 
phase transitions induced by the FFLO phase in Pauli limited multi-band iron-based 
superconductors at low temperature and high magnetic field. 

2. Model and theoretical background 

IBSC are chemical compounds possessing layered structure, with a characteristic Fermi 
surface (FS). In order to explain the FS features in IBSC, various two-orbital [39], 
three-orbital [40,41] and five-orbital [42,43] tight binding models have been proposed. 
The abundance of available models is induced by the relatively complicated band 
structure of IBSC materials, which is strongly dependent on the chemical doping [44]. 
To realistically describe e.g. FeAs layers, five 3d-orbitals of iron ions need to 
be retained in the model. However, the band structure calculations suggest the 
importance of itinerant electrons from the d xz and d yz orbitals. Influential in the 
formation of the FS are also the d xy orbitals. For this reason, in our calculations we 
use for simplicity the three-band model proposed by M. Daghofer et al. [40,41]. 

In general, the momentum-dependent tight-binding non-interacting Hamiltonian 
of the multi-orbital IBSC in orbital space is be given by: 

tfo = ee(^ - (p + ah)S a ^j cl krT cp ka , (1) 

a/3 ka 

where c a k a (c[, ) annihilates (creates) an electron with momentum k and spin er in 
the orbital a. The hopping matrix elements correspond to the kinetic energy 
of a particle with momentum k changing the orbital from fd to a, they are given by 
the effective tight-binging model of the two-dimensional FeAs planes in the selected 
model. We use the model of IBSC proposed by M. Daghofer et al. in Ref. [40] and 
improved in Ref. [41]. Beyond the d xz and d yz orbitals, the model also accounts for 
the d xy orbital: 

T^ 1 = 2f 2 cos k x + 2t\ cos k y + 4f 3 cos k x cos k y 

+ 2tu(cos(2fc x ) — cos(2 k v )) (2) 

+ 4f 12 cos(2 k x ) cos(2 k y ), 

T| 2 = 2f i cos k x + 21 2 cos k y + 41 3 cos k x cos k y 

— 2f 11 (cos(2fe r ) — cos(2 k y )) (3) 

+ 4f 12 COS^fca,) cos(2 ky), 

T£ 3 = eo + 2t 5 (cos k x + cos k y ) + 4te cos k x cos k y 

+ 2t 9 {cos(2k x ) + cos(2 k y )) (4) 

+ 4fi 0 (cos(2fc x ) cos k y + cos k x cos(2 k y )), 

= T/p = 4ti sin k x sin k y , 


Tl 2 


( 5 ) 
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T^ 3 = Tfe 1 = 2it 7 sin k x + 4 its sin k x cos k y , (6) 

T| 3 = f| 2 = 2it 7 sin k v + 4its sin k y cos k x . (7) 

In Ref. [41] the hopping parameters are given in electron volts as: = —0.08, 


t 2 = 0.1825, t 3 = 0.08375, t 4 = -0.03, t 5 = 0.15, t 6 = 0.15, t 7 = -0.12, t 8 = 0.06, 
tg = 0.0, tig = —0.024, t\i = —0.01, tig = 0.0275 and eo = 0.75. The average number 
of particles in the system n = 4 is attained for p = 0.4748. 

The band structure of the IBSC model can be reconstructed from the kinetic tight- 
binding Hamiltonian in the orbital representation Hg via the unitary transformation 
H' 0 = WHgU [40]. Then Hg = J2 s ka E Ekc rd\ ktT d EkcT . Here d Ek(J ( d \ ktT ) annihilates 
(creates) an electron in band e. The total number of particles in the system 
J2akc r c ^aka c aka = J2 e ka d\ ka d kea is adjusted by the chemical potential p. We 
neglected orbital effects, which is equivalent to assuming that the external magnetic 
field h is parallel to the FeAs layers. 

For simplicity and readability we assume the existence of only the superconducting 
phase in the system. However, it should be had in mind that in many IBSC 
systems superconductivity can coexist with magnetic order. [45] In this paper, without 
specifying the mechanisms responsible for the forming of superconducting phases, we 
introduce a superconducting pairing between the gwasi-particles in bands e, which 
is a good approximation in the limit of weak or vanishing inter-band pairing [46]. 
Superconducting states with non-zero total momentum of Cooper pairs (TMCP) can 
be described by the phenomenological effective Hamiltonian: 

H'sc = + H c ■) > ( 8 ) 

ke 

where A ek = A E rj(k) = U e p(k)(d e -k+q e id e k 7 ) is the SOP in band e for the TMCP q E 
and amplitude A E . The structure factor given by r)(k) captures the symmetry of the 
SOP, related to the effective interaction in real space [26,47]. H' sc in band space is 
the reformulation of the interacting Hamiltonian in orbital space [40] Similarly to the 
two-band model, the SOP in the band representation can be transformed to orbital 
one. [48] Moreover, the interband SOPs with different values of A ek in every band e 
correspond to the existence of the intra- and interorbital SOPs in the system. 

The total Hamiltonian H = H' 0 + H' sc , formally describes a system with a three 
independent bands. For this reason, the eigenvalues of H in the band representation 
are given by standard Bogoliubov transformation [49,50]: 

= n ± ^«) 2 + |A efc P Vfc, (9) 

where dj = {E ek -f ± E Et _ k+qe ±) /2. The grand canonical potential can be calculated 
explicitly from its definition Q, = A:sTln{Tr[exp(— H/kgT)]}, which for given (fixed) 
parameter h and T can be treated as function of the SOP A e and TMCP q E in each 
bands. 

2.1. The entropy and specific heat calculation 

In case of intraband superconductivity, the grand canonical potential is given as 
Cl = where: 


( 10 ) 
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|A efc | s 

U F 


Detail calculation can be found in Ref. [26, 48-50]. From the thermodynamical 
potential, we can determine the entropy S = —d£l/dT and superconducting specific 
heat at temperature T as C = — Td 2 £l/dT 2 , where the is the grand canonical 
potential. The entropy [51] is S = —dfl/dT, where: 


dT 



dflA (dtl e \ dA e 
dTj e + {d Aj e dT 



dq e ' 

V dq e , 

' e ST. 


( 11 ) 


where subscript e labels the equilibrium values of the SOPs A e and TMCPs q e . From 
the equilibrium condition, we have: 9fi e /5A e | e = d£l e /dq e \e = 0 for all e, since 
f2(A e , q £ ) is at a minimum. Hence S = — dtt s /dT\e, or a priori: 


5 = E 


eh 


\ot 

4f/( a: 


eh) 


) + Ub In ( 1 + exp ( - 


A' 


eh 


k B T 


( 12 ) 


where f(E) is the Fermi-Dirac distribution. The specific heat is then defined in the 
usual manner by: C = TdS/dT\ n ^,h,v = — T9 2 fl/<9T 2 | e . It should be noted that 
the grand potential (and also SOPs and TMCPs) depend on temperature in a non¬ 
trivial manner, which forces the calculation of S and C to be carried out by numerical 
derivatives. 



Figure 1. (Color on-line) (a-c) Superconducting intra-band order parameter |A e | 
(energetic gap amplitude) for each band e, in function of the external magnetic 
field for different temperatures: T equal 0.10 (red dotted line), 0.08 Tq (blue 
solid line), and 0.06 (green dashed line), gray solid line denotes results for 
T~0. (d) Cooper-pair total momentum q £ in each band e, as function of the 
applied magnetic field h. The lines end at the critical field of the corresponding 
band. The temperature dependence is negligible in the considered regime. Pink 
points mark the phase transitions for the different temperatures (given above). 



































Multiple phase transitions in Pauli limited iron-based superconductors 6 

3. Numerical results and discussion 

The global ground state can be obtained, at fixed values of the parameters 
(temperature T and magnetic field h), from the minimization of the grand canonical 
potential with respect to the SOP A e and TMCP q £ , at the same time determining 
the optimal values of the latter parameters. All calculations have been performed on 
NVIDIA GPUs, in momentum space on a square lattice grid k x x k y = 10000 x 10000, 
using the algorithm described in Ref. [50]. The following are predictions for s-wave 
symmetry of the SOP, however other symmetries generate analogous results. 

We assumed a different effective attractive intraband pairing in each band. It is 
found by seeking the simultaneous disappearance of the superconducting BCS phase 
at critical temperature Tc, equal to 5 meV (~55 K). For this set of parameters, the 
BCS critical magnetic h® cs is also uniform in the bands, with a value of about 6 
meV (~103 T). Although these relatively large values can be found in some class of 
IBSC [52], we focus on generic features of IBSC in LTHM regime. 

Above h .£ CS , at temperatures below some characteristic T + , the FFLO phase 
arises [6]. In this phase the SOPs decrease in increasing external magnetic field 
(Fig. 1 a-c). Moreover, in general the TMCP depend on the size of the splitting 
between the Fermi surfaces for electrons with spin up and down, the source of which 
is the Zeeman effect. Raising the external magnetic field increases the splitting, which 
leads to greater TMCP (Fig. l.d) [48], which is also true in other systems [53]. The 
TMCP dependence on the temperature is however weak. 

In the case of BCS, we can find a typical h T phase diagram (Fig. 2. a), in which 
superconductivity in every band disappears at the same external magnetic field [48]. 
Conversely for the FFLO phase, a different critical magnetic field h^ LO [26, 50] 
determines the superconducting behavior in each band, as shown in Fig. 2. a. Because 
of this, calorimetric experiments in LTHM regime should display multiple phase 
transitions , in the form of a group of peaks in specific heat, shown as arrows in Fig. 2.b. 
The first group (or extended peak) can be connected with the transition from the BCS 
to the FFLO phase (first arrow from the left in Fig. 2.b), while the second group arises 
due to transitions inside the FFLO phase (blue and green line in Fig. 2. a), and to the 
final transition from superconducting to the normal state (red line in Fig. 2. a). To 
conclude, the measured total specific heat is characterized by a multiply discontinuous 
shape, the black line in Fig. 2.b. It should be noted, that similar effect can be found in 
systems with relatively significant finite-size effects. [54-56] However, multiple phase 
transitions in this case (and the relative discontinuities in the specific heat), are due to 
abrupt changes in the TMCP q e , and not the result of the disappearance of the FFLO 
phase in consecutive bands. This same effect can be expected in two dimensional 
square lattice, where a growing magnetic field favors FFLO phases with a greater 
number of inequivalent momenta entering the TMCP [57]. 

Usually the superconductivity in IBSC is described by an s±-wave symmetry, 
where the gap changes its sign between the hole and electron pockets of the Fermi 
surface [58,59]. However, d-wave symmetry can be also observed [60,61]. Hole doping 
of IBSC can lead to the transition of the gap symmetry from s±-wave to d-wave [62-64]. 
Our main result is not affected, with similar results for other symmetries than s-wave - 
the FFLO is the ground state in the LTHM regime in every band. However, a different 
SOP symmetry could influence the shape of specific heat in function of temperature 
(Fig. 2). Moreover, because the band structure of the IBSC strongly depends on 
doping, more realistic results require specific models to the chemical compound. 
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For different compounds, the quantity of specific heat peaks can depend on the 
number of bands forming the Fermi surface and supporting the FFLO state in the 
LTHM regime, with the detailed band structure influencing only quantitatively the 
jump heights. 

4. Summary and final remarks 

In this paper we show that the presence of the FFLO phase, in multi-band materials 
like IBSC, can be experimentally ascertained by the appearance of multiple phase 
transitions , which is in turn manifested by multi-discontinuities in the shape of the 
specific heat. In the context of Pauli limited superconductors, we can speak about 
two scales of temperature. The former is the critical temperature T c, at which 



0.0 0.4 0.8 1.2 

T[meV] 


Figure 2. (Color on-line) (a) Magnetic field h - temperature T phase diagram. 
Solid lilies denote the phase transitions between BCS state and normal or FFLO 
state, while red, green and blue lines phase transitions inside the FFLO state, 
described in the main text, (b) Specific heat along the cut line (violet dotted 
line in panel a). Bold black line is the total specific heat while red, green and 
blue solid lines are the partial specific heats for the first, second, and third band 
respectively. Dotted lines display the form of the specific heat for the normal 
state. Arrows mark the locations of phase transitions. 
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the superconductivity vanishes in zero external magnetic field. The latter is the 
temperature T + , in which phase transitions from the superconducting to normal state 
change kind, from first to second order. is reported as 2.3 K and 3.5 K for 
CeCoIns and KFe 2 As 2 respectively, while in both cases T + can be approximated 
as 0.31 T c [8, 19]. In CeCoIn 5 the specific heat displays an additional anomaly 
within the superconducting state at a temperature ~ 300 rnK (~ 0.12 T<y) [7]. These 
experimental results are interpreted as evidence for the existence of the FFLO phase. 
Similar behavior is observed in organic superconductors [14]. 

Relevant experimental data have been presented in the literature for KFe 2 As 2 [65, 
66], but not at sufficiently low temperatures. The iron-based KFe 2 As 2 Pauli-limited 
superconductors have critical magnetic field near 5 T [19]. However also for fields above 
this value, an anomalous shape for the specific heat can be observed in LTHM regime 
(for temperatures below 1 K) [65,66]. The anomalies have been related to a meta- 
magnetic transition in high-field heavy-fermion Celrlns [67] , which is explained by the 
magnetization process in the antiferromagnet with helical spin structure. However, 
the heavy-fermion CeCoIns additional phase transition in the LTHM regime inside the 
superconducting phase can be the manifestation of the emergence of incommensurate 
spin density wave (SDW) [68]. The latter mechanism is strongly connected with 
the FFLO phase, because an existing SDW increases the tendency of the system 
toward the creation of an FFLO phase [69]. To verify this hypothesis, accurate 
calorimetric measurements are required in LTHM regime. This phenomenon requires 
further experimental and theoretical studies. 
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